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Abstract: The highly unsaturated binuclear butadiene iron carbonyls (C4H6)2Fe2(CO)^ 
(fz = 2, 1) have been examined using density functional theory. For (C4H6)2Fe2(CO)^ (n = 2, 1), 
both coaxial and perpendicular structures are found. The global minima of 
(C4H6)2Fe2(CO)^ (n = 2, 1) are the perpendicular structures 2Q-1 and lQ-1, respectively, 
with 17- and 15-electron configurations for the iron atoms leading to quintet spin states. 
The Fe=Fe distance of 2.361 A (M06-L) in the (C4H6)2Fe2(CO)2 structure 2Q-1 suggests a 
formal double bond. The Fe=Fe bond distance in the (C4H6)2Fe2(CO) structure lQ-1 is 
even shorter at 2.273 A (M06-L), suggesting a triple bond. Higher energy (C4H6)2Fe2(CO)^ 
(n = 2, 1) structures include structures in which a bridging butadiene ligand is bonded to 
one of the iron atoms as a tetrahapto ligand and to the other iron atom through two agostic 
hydrogen atoms from the end CH2 groups. Singlet (C4H6)2Fe2(CO) structures with formal 
Fe-Fe quadruple bonds of lengths --2.05 A were also found but at very high energies 
(-47 kcal/mol) relative to the global minimum. 

Keywords: iron carbonyls; iron-iron bonding; agostic hydrogen atom; 
metal-olefin complexes 
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1. Introduction 

The chemistry of metal carbonyl complexes of acyclic hydrocarbons dates back to the 1930 
discovery by Reihlen et al. [1] of the mononuclear butadiene iron tricarbonyl complex, C4H6Fe(CO)3 
by the reaction of butadiene with iron pentacarbonyl at elevated temperatures. The proposed tetrahapto 
bonding of the butadiene ligand to the Fe(CO)3 unit in this complex was confirmed in 1963 by Mills 
and Robinson [2] using X-ray crystallography at -40 °C (Figure lA). In addition, in 1962 Murdoch 
and Weiss [3] used the reaction of butadiene with Fe2(CO)9 at room temperature to synthesize the 
tetracarbonyl (r| -C4H6)Fe(CO)4 in which only one of the two C=C double bonds of the butadiene 
ligand is bonded to the iron atom. An additional product from the latter reaction was the binuclear 
complex C4H6[Fe(CO)4]2 in which each C=C double bond of the butadiene ligand is bonded to a 
separate Fe(C0)4 unit with the iron atoms much too far apart for any kind of direct iron-iron bond. 
However, no binuclear butadiene iron carbonyl derivatives with short iron-iron distances suggesting 
iron-iron bonds have been synthesized. In order to assess the possibilities for binuclear iron carbonyl 
derivatives with iron-iron bonds we have performed a density functional theory (DFT) study on 
possible structures for (r|'^-C4H6)2Fe2(CO)„ {n = 5, 4, 3), predicted to have structures with formal Fe-Fe 
single bonds, Fe=Fe double bonds, and Fe=Fe triple bonds, respectively [4]. In general, the lowest 
energy structures for these (r|'^-C4H6)2Fe2(CO)„ derivatives were found to be coaxial structures in 
which each metal atom is bonded to a single butadiene ligand (Figure IB). 

Figure 1. Structure of (r|'^-C4H6)Fe(CO)3 (A) and two general structure types for binuclear 
metal carbonyls, coaxial (B) and perpendicular (C). 




A B C 



This paper reports a DFT study of the still more highly unsaturated (C4H6)2Fe2(CO)„ (n = 2, 1) 
derivatives. Such systems are interesting since strict adherence to the 18 -electron rule suggests that 
these highly unsaturated derivatives might provide examples of very short formal iron-iron quadruple 
and quintuple bonds. However, for such systems containing two or fewer carbonyl groups, alternative 
perpendicular structures are possible in which the butadiene ligands bridge the pair of iron atoms 
(Figure IC). Structures of both types were found in this work. Quintet spin state structures were found 
to be the lowest energy structures for both (C4H6)2Fe2(CO)2 and (C4H6)2Fe2(CO). In addition, 
interesting structures were found with agostic hydrogen atoms bridging Fe-C bonds. 
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2. Theoretical Methods 

Density functional theory (DFT) including electron correlation effects has been shown to be a 
powerful and effective computational tool in organotransition metal chemistry [5-19]. Three DFT 
methods were used for this work. The first functional was BP86, which combines Becke's 1988 
exchange functional (B) with Perdew's 1986 gradient corrected correlation functional method (P86), 
and usually provides better vibrational frequencies [20,21]. The second DFT method was B3LYP, 
which is the hybrid HF/DFT functional using the combination of the three -parameter Becke functional 
(B3) [22] with the Lee- Yang-Parr (LYP) generalized gradient correlation functional [23]. The third 
was a hybrid meta-GGA DFT method, M06-L, developed by Truhlar's group [24]. Recently Truhlar's 
group made much progress to the development of improved exchange-correlation functionals that are 
essential for expanding the applicability of Kohn-Sham DFT, such as the M06 suite. Thus M06-L was 
constructed using three strategies: constraint satisfaction, modeling the exchange-correlation hole, and 
empirical fits. They concluded that M06-L is one of the best functionals for the study of 
organometallic and inorganic thermochemistry, and is the best functional for transition metal 
energetics. In comparing the first two DFT methods Reiher and collaborators found that B3LYP 
always overestimates the energy of high-spin states and BP86 overestimates the energies of low-spin 
states for a series of the Fe(II)-S complexes [25]. In the present study, we found that the M06-L 
method predicts an intermediate energy difference, anticipated to be closer to the experimental. We 
therefore adopt the energy order predicted by the M06-L method, but list the BP86 and B3LYP results 
in the Supporting Information. 

Basis sets have been chosen to provide continuity with a body of existing research on 
organometallic compounds. Fortunately, DFT methods are far less basis set sensitive than higher-level 
methods such as coupled cluster theory. In this work, the double-^ plus polarization (DZP) basis sets 
used for C and O add one set of pure spherical harmonic d functions with orbital exponents 
ad(C) = 0.75 and ad(0) = 0.85 to the Huzinaga Dunning standard contracted DZ sets and are denoted 
as (9s5pld/4s2pld) [26,27]. For H, a set of p polarization functions ap(H) = 0.75 is added to the 
Huzinaga Dunning DZ sets. For Fe, in our loosely contracted DZP basis set, the Wachters' primitive 
set [28] is used after being augmented by two sets of p functions and one set of d functions and then 
contracted using the method of Hood, Pitzer, and Schaefer [29]. This basis set is denoted as 
(14sllp6d/10s8p3d). 

The geometries of the structures were fully optimized using the Gaussian09 program [30] with the 
three selected DFT methods and with the indicated DZP basis set. The vibrational frequencies were 
determined by evaluating analytically the second derivatives of the energy with respect to the nuclear 
coordinates at the same levels. The corresponding infrared intensities were also evaluated analytically. 
The fine grid (75, 302) was the default for evaluating integrals numerically, and the tight designation 
was the default for the energy convergence, as well as the tight option for the geometry optimizations. 
In some cases, the finer grid (120, 974) was used for investigating small imaginary vibrational 
frequencies. Natural bond orbital (NBO) analyses [31] used the DZP BP86 method with the NBO 3.1 
version attached in the Gaussian03 program. 

The optimized (C4H6)2Fe2(CO)^ {n = 2, 1) structures are depicted in Figures 2 to 7. In these figures, 
the upper and lower distances were obtained by the M06-L and BP86 method, respectively. The 
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structures are designated as nX-m, where n stands for the number of CO groups, X designates the spin 
state using S for singlets, T for triplets and Q for quintets, and m orders the structures according to 
their relative energies. Note that, although the singlets, the triplets and the quintets are discussed in 
separate sections, the relative energies are considered together on the basis of the number of carbonyls. 
The M06-L method appears to predict the better singlet-triplet splittings. 

The relative energies corrected for zero-point energies are listed in the Supporting Information, 
where computed enthalpies and free energies are also given. These relative free energies agree within 
2 kcal/mol with the relative electronic energies. 

3. Results and Discussion 

3.1. (C4H6)2Fe2(CO)2 

Three types of stationary points, namely, coaxial structures, perpendicular structures, and deformed 
coaxial structures, have been found for (C4H6)2Fe2(CO)2. The global minimum of the (C4H6)2Fe2(CO)2 
is the quintet structure 2Q-1 according to the relative energies listed in Tables 1, 2, and 3. 

3.1.1. Quintet (C4H6)2Fe2(CO)2 Structures 

Three quintet structures were found for (C4H6)2Fe2(CO)2 (Figure 2 and Table 1). The perpendicular 
structure 2Q-1 with a bridging CO group and two bridging butadiene ligands is the global minimum. 

o 

The Fe=Fe distance in 2Q-1 is predicted by M06-L to be 2.361 A. This can be interpreted as a formal 
Fe=Fe double bond to give one iron atom a 17-electron configuration and the other iron atom a 
15-electron configuration. This can correspond to a quintet spin state. 

Figure 2. Quintet structures for (C4H6)2Fe2(CO)2. The upper bond distances were obtained 
by the M06-L method and the lower bond distances by the BP86 method. 




2Q-1,C, 



2Q-2, Ci 



2Q-3, Czh 
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Table 1. Fe-Fe distances (A), HOMO-LUMO energies (E, in hartree), HOMO-LUMO 
gaps (in eV), total energies (E, in hartree), relative energies (AE, in kcal/mol, numbers of 
imaginary frequencies (Nimag) and spin expectation values {S ) for the quintet 
(C4H6)2Fe2(CO)2 structures with the M06-L method. 



2Q-1 (C) 



2Q-2 (CO 



2Q-3 (C2h) 



M06-L Fe-Fe 



2.361 
-0.17183 
-0.09384 



2.311 
-0.17530 
-0.10290 

1.97 
3066.14371 
1.5 



2.342 
-0.17208 
-0.11149 

1.65 
3066.13340 
8.0 



HOMO(a) 

LUMO(a) 

gap/eV 

E 

AE 



2.12 
3066.14608 
0.0 



Nimag 



none 



6.30 



none 



6.31 



none 



6.63 



A deformed coaxial (C4H6)2Fe2(CO)2 structure 2Q-2 with one terminal CO group, one semibridging 
carbonyl group, and one bridging butadiene ligand is predicted to lie only 1.5 kcal/mol (M06-L) in 
energy above the global minimum 2Q-1 (Figure 2 and Table 1). The bridging butadiene ligand in 2Q-2 
is bonded as a trihapto ligand to one iron atom and as a monohapto ligand to the other iron atom. The 
Fe=Fe bond distance in 2Q-2 of 2.311 A (M06-L) is close to that in 2Q-1 and can likewise be 
interpreted as a formal double bond. Again this gives one iron atom a 15 -electron configuration and 
the other iron atom a 17-electron configuration, which can correspond to a quintet spin state. 

The coaxial (C4H6)2Fe2(CO)2 structure 2Q-3 has two bridging carbonyl groups and lies 8.0 kcal/mol 
in energy above 2Q-1. The Fe=Fe bond distance in 2Q-3 is 2.342 A (M06-L), which is similar to those 
in 2Q-1 and 2Q-2 and thus can correspond to a formal double bond. This gives each iron atom in 2Q-3 
a 16-electron configuration, which can correspond to a quintet spin state. 

3.1.2. Triplet (C4H6)2Fe2(CO)2 Structures 

Three triplet structures were found for (C4H6)2Fe2(CO)2 (Figure 3 and Table 2). The lowest energy 
of these triplet structures, namely 2T-1 lying 6.8 kcal/mol in energy above the quintet global minimum 
2Q-1, is a coaxial structure with two bridging carbonyl groups and two terminal butadiene ligands. 
The Fe^Fe bond distance of 2.209 A (M06-L) in 2T-1 is --0.1 A shorter than that in the similar quintet 
spin state structure 2Q-3 and thus can be interpreted as a formal triple bond. This gives each iron atom 
in 2T-1 a 17-electron configuration consistent with a binuclear triplet. 

The next triplet (C4H6)2Fe2(CO)2 structure, namely 2T-2, lying 8.4 kcal/mol above the quintet 
global minimum 2Q-1, is a deformed coaxial structure with one bridging butadiene ligand, one 
terminal butadiene ligand, and two terminal carbonyl groups (Figure 3 and Table 2). The bridging 
butadiene ligand in 2T-2 is bonded to one of the iron atoms as a tetrahapto ligand. In addition, the two 
terminal hydrogen atoms of this butadiene ligand are agostic hydrogen atoms bonding to the other iron 

o o 

atom through C-H-Fe bridging units with Fe-C distances of -2.2 A and Fe-H distances of -1.9 A. 
These C-H-Fe units are predicted to exhibit abnormally low v(C-H) frequencies of 2515 and 
2561 cm~\ The Fe-Fe distance in 2T-2 is relatively long at 2.433 A (M06-L) and can thus be 
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considered as a formal single bond. This gives each iron atom in 2T-2 the 17-electron configuration 
for a binuclear triplet. 



Figure 3. Triplet structures for (C4H6)2Fe2(CO)2. The upper bond distances were obtained 
by the M06-L method and the lower bond distances by the BP86 method. 




2T.1, C2 2T.2, Ct 2T-3, C2 



Table 2. Fe-Fe distances (A), HOMO-LUMO energies (E, in hartree), HOMO-LUMO 
gaps (in eV), total energies (E, in hartree), relative energies (AE, in kcal/mol), numbers of 
imaginary frequencies (Nimag) and spin expectation values (5 ) for the triplet 
(C4H6)2Fe2(CO)2 structures with the M06-L method. 





2T-1 (C2) 


2T-2 (Ci) 


2T-3 (C2) 


M06-L Fe-Fe 


2.209 


2.433 


2.295 


HOMO(a) 


-0.19453 


-0.17557 


-0.15910 


LUMO(a) 


-0.12453 


-0.10739 


-0.10465 


gap/eV 


1.90 


1.86 


1.48 


E 


-3066.13525 


-3066.13271 


-3066.1265 


AE 


6.8 


8.4 


12.3 


Nimag 


none 


none 


none 




2.20 


2.12 


2.21 



The perpendicular triplet (C4H6)2Fe2(CO)2 structure 2T-3 has two bridging butadiene ligands and 
two terminal carbonyl groups and lies 12.3 kcal/mol (M06-L) in energy above the global minimum 
2Q-1 (Figure 3 and Table 2). The Fe=Fe distance of 2.295 A (M06-L) in 2T-3 can correspond to a 
formal triple bond to give each iron atom a 17-electron configuration for a binuclear triplet. 



3.1.3. Singlet (C4H6)2Fe2(CO)2 Structures 

Three distinct singlet structures are obtained for (C4H6)2Fe2(CO)2. However, they are all high 
energy structures lying from 19 to 30 kcal/mol above the 2Q-1 global minimum (Figure 4 and Table 3). 
The lowest energy singlet (C4H6)2Fe2(CO)2 structure 2S-1, lying 19.1 kcal/mol above 2Q-1, is very 
similar to the triplet structure 2T-2. Thus structure 2S-1 has a bridging butadiene ligand connected to 
one iron atom as a tetrahapto ligand and to the other atom through two non-equivalent C-H-Fe bridges 
with Fe-C distances of 2.166 and 2.431 A and Fe-H distances of 1.844 and 2.113 A (M06-L). The 
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strikingly lower v(C-H) vibrational frequencies 2392 cm"^ and 2912 cm"^ (BP86) and the longer C-H 
bond distances involving these agostic hydrogens confirm the weaker C-H bonds. One of the carbonyl 
groups in 2S-1 is a nearly symmetrical bridging carbonyl group with Fe-C distances of 1.878 and 

o 

1.948 A (M06-L). The remaining carbonyl group and butadiene ligand in 2S-1 are both terminal 
groups. The predicted Fe=Fe bond length of 2.344 A (M06-L) in 2S-1 is --0.1 A longer than the formal 
Fe-Fe single bond in 2T-2 and thus can be considered to be a formal double bond, thereby giving both 
iron atoms the favored 18-electron configuration in 2S-1. 

The C2h coaxial singlet (C4H6)2Fe2(CO)2 structure 2S-2 with two bridging carbonyl groups and 
terminal butadiene ligands is predicted to lie 23.9 kcal/mol (M06-L) in energy above the global 

o 

minimum 2Q-1 (Figure 4 and Table 3). The Fe=Fe bond distance of 2.327 A (M06-L) in 2S-2 is 
similar to that in 2S-1 and thus can be assigned to a formal double bond. This gives each iron atom in 
2S-2 a 16-electron configuration, which can relate to a binuclear singlet. 

Figure 4. Singlet structures for (C4H6)2Fe2(CO)2. The upper bond distances were obtained 
by the M06-L method and the lower bond distances by the BP86 method. 




The C2 singlet (C4H6)2Fe2(CO)2 structure 2S-3 is a perpendicular structure lying 30.0 kcal/mol 
(M06-L) above the 2Q-1 global minimum (Figure 4 and Table 3). Both butadiene ligands in 2S-3 are 
bridging ligands and both carbonyl ligands in 2S-3 are terminal ligands. The Fe=Fe distance of 
2.325 A in 2S-3 is similar to those in 2S-1 and 2S-2 and likewise can correspond to a formal double 
bond. This gives each iron atom in 2S-3 a 16-electron configuration suggesting a vacant coordination 
site on each iron atom. This is consistent with the geometry of the carbonyl groups in 2S-3. 

The overall energy order for (C4H6)2Fe2(CO)2 structures investigated in this section is 2Q-1 < 2Q-2 
< 2T-1 < 2Q-3 - 2T-2 < 2T-3 < 2S-1 < 2S-2 < 2S-3 by M06-L. This suggests that (C4H6)2Fe2(CO)2 
prefers high spin state structures. The energy gaps between the HOMOs and LUMOs in Tables 1, 2 
and 3 for 2S-1 2T-2 2Q-2 are 1.01, 1.86 and 1.97 eV, respectively, which increase 
monotonically with the increase of the spin multiplicity. The same trend is also found for 2S-2 2T-1 

2Q-3 and 2S-3 2T-3 2Q-1. 
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Table 3. Fe-Fe distances (A), HOMO-LUMO energies {E, in hartree), HOMO-LUMO 
gaps (in eV), total energies {E, in hartree), relative energies (A£', in kcal/mol) and numbers 
of imaginary frequencies (Nimag) for the singlet (C4H6)2Fe2(CO)2 structures with the 



M06-L method. 




2S-1 (CO 


2S-2 (C2h) 


2S-3 (C2) 


M06-L Fe-Fe 


2.344 


2.327 


2.325 


HOMO 


-0.15788 


-0.17415 


-0.13105 


LUMO 


-0.12068 


-0.13219 


-0.11009 


gap/eV 


1.01 


1.14 


0.57 


E 


-3066.11570 


-3066.10797 


-3066.09822 


AE 


19.1 


23.9 


30.0 


Nimag 


none 


none 


none 



3.2. (C4H6)2Fe2(CO) 

Two types of stationary points for (C4H6)2Fe2(CO), namely, coaxial and perpendicular structures, 
are shown in Figures 5, 6 and 7. The global minimum of the (C4H6)2Fe2(CO) is the quintet structure 
lQ-1 according to the relative energies listed in Tables 4, 5, and 6. 

3.2.1. Quintet (C4H6)2Fe2(CO) Structures 

Two quintet structures were found for (C4H6)2Fe2(CO) (Figure 5 and Table 4). The perpendicular 
(C4H6)2Fe2(CO) global minimum lQ-1 has two bridging butadiene ligands and a terminal carbonyl 
group. The Fe=Fe distance of 2.273 A (M06-L) can correspond to a formal triple bond thereby giving 
the iron atom bearing the carbonyl group a 17-electron configuration but the other iron atom only a 
15-electron configuration. This is consistent with a quintet spin multiplicity. 

Figure 5. Quintet structures for (C4H6)2Fe2(CO). The upper bond distances were obtained 
by the M06-L method and the lower bond distances by the BP86 method. 




The other quintet (C4H6)2Fe2(CO) structure lQ-2 has a bridging CO group and lies 12.9 kcal/mol 
above the global minimum lQ-1 (Figure 5 and Table 4). The Fe=Fe distance of 2.283 A (M06-L) in 
lQ-2 is similar to that in lQ-1 and likewise can correspond to a formal triple bond. This gives both 
iron atoms a 16-electron configuration, which can correspond to a binuclear quintet spin state. 
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Table 4. Fe-Fe distances (A), HOMO-LUMO energies {E, in hartree), HOMO-LUMO 



gaps (in eV), total energies {E, in hartree), relative energies (A£', in kcal/mol, numbers of 
imaginary frequencies (Nimag) and spin expectation values (5 ) for the quintet 
(C4H6)2Fe2(CO) structures with the M06-L method. 



iQ-1 (Ci) 



lQ-2 (Ci) 



M06-L 



Fe-Fe 

HOMO(a) 

LUMO(a) 

gap/eV 

E 

AE 



2.290 
-0.16350 
-0.07858 

2.31 
■2952.77563 
0.0 



2.283 
-0.17107 
-0.10795 



1.72 
2952.75505 
12.9 



Nimag 



none 



6.32 



none 



6.40 



3.2.2. Triplet (C4H6)2Fe2(CO) Structures 

Two triplet (C4H6)2Fe2(CO) structures were found to lie 12 to 24 kcal/mol in energy above the 
quintet global minimum lQ-1 using the M06-L method (Figure 6 and Table 5). The perpendicular 
triplet (C4H6)2Fe2(CO) structure lT-1, lying 11.6 kcal/mol above lQ-1, has two bridging butadiene 
ligands and a terminal carbonyl group. Considerable spin contamination was found for lT-1. Thus the 
spin expectation value (S^) = 2.82 for lT-1 as compared with the ideal 2.0. Indeed, a more stable 
quintet (C4H6)2Fe2(CO) structure lQ-1 is found with an {S ) value within 10% of the ideal 6.0. The 

o 

Fe=Fe distance of 2.291 A (M06-L) is very similar to the Fe=Fe triple bond distances in the quintet 
structures lQ-1 and lQ-2 and thus likewise can correspond to a formal triple bond. This gives the iron 
atom in lT-1 bearing the carbonyl group a 17-electron configuration but the other iron atom only a 
15-electron configuration. This can correspond to a binuclear triplet with a vacant coordination 
position on the iron atom with only a 15-electron configuration. The Cs coaxial triplet (C4H6)2Fe2(CO) 
structure lT-2, lying 24.2 kcal/mol (M06-L) above lQ-1 has a bridging carbonyl group but terminal 

o 

butadiene ligands. The predicted Fe=Fe distance of 2.200 A (M06-L) is similar to the Fe=Fe distances 
in the other quintet and triplet (C4H6)2Fe2(CO) structures and likewise can correspond to a formal 
triple bond. This gives each iron atom in lT-2 a 16-electron configuration. 

Figure 6. Triplet structures for (C4H6)2Fe2(CO). The upper bond distances were obtained 
by the M06-L method and the lower bond distances by the BP86 method. 




1T-1,C 



's 



1T-2, Cs 
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3.2.3. Singlet (C4H6)2Fe2(CO) Structures 

Two singlet low-lying (C4H6)2Fe2(CO) structures were found (Figure 7 and Table 6) but at very 
high energies relative to the corresponding quintet and triplet structures. Thus the Cs perpendicular 
singlet structure lS-1 lies 46.7 kcal/mol above the quintet global minimum lQ-1. Structure lS-1 has a 
terminal carbonyl group and bridging butadiene ligands. The Fe-Fe distance in lS-1 of 2.051 A is 
--0.2 A shorter than the Fe=Fe distances in the quintet and triplet (C4H6)2Fe2(CO) structures and thus 
can correspond to a formal quadruple bond. This gives the iron atom in lS-1 bearing the carbonyl 
group the favored 18-electron configuration but the other iron atom only a 16-electron configuration. 

Table 5. Fe-Fe distances (A), HOMO-LUMO energies (£, in hartree), HOMO-LUMO 
gaps (in eV), total energies {E, in hartree), relative energies (AE", in kcal/mol, numbers of 
imaginary frequencies (Nimag) and spin expectation values for the triplet 

(C4H6)2Fe2(CO) structures using the M06-L method. 







lT-1 (Q 


lT-2 (Q 


M06-L 


Fe-Fe 


2.291 


2.200 




HOMO(a) 


-0.15400 


-0.164149 




LUMO(a) 


-0.102794 


-0.093119 




gap/eV 


1.39 


1.93 




E 


-2952.75716 


-2952.73705 




A£ 


11.6 


24.2 




Nimag 


none 


none 




(5^) 


2.82 


2.22 



Figure 7. Singlet structures for (C4H6)2Fe2(CO). The upper bond distances were obtained 
by the M06-L method and the lower bond distances by the BP86 method. 




1S-1,Cs 1S-2, Cs 



The Cs singlet coaxial (C4H6)2Fe2(CO) structure lS-2 is also a high energy structure, lying 
47.4 kcal/mol above the global minimum lS-1 (Figure 7 and Table 6). Structure lS-2 has terminal 
butadiene ligands and a bridging carbonyl group. The very short Fe-Fe distance of 2.039 A (M06-L) in 
lS-2 is similar to that in lS-1 and thus likewise can correspond to a formal quadruple bond. This gives 
one iron atom in lS-2 the favored 18-electron configuration but the other iron atom only a 16-electron 
configuration. This asymmetry in the electron count on the iron atoms in lS-2 is reflected in a different 
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arrangement of the rj'^-butadiene ligands on each iron atom. Thus, the "right" iron atom in lS-2 as 
depicted in Figure 7 appears to have a vacant coordination site and thus can correspond to the iron 
atom with a 16-electron configuration. 

Table 6. Fe-Fe distances (A), HOMO-LUMO energies {E, in hartree), HOMO-LUMO 
gaps (in eV), total energies (£", in hartree), relative energies (AE", in kcal/mol, and numbers 
of imaginary frequencies (Nimag) for the singlet (C4H6)2Fe2(CO) structures with the 
M06-L method. 







lS-1 (Q 


lS-2 (Q 


M06-L 


Fe-Fe 


2.051 


2.039 




HOMO 


-0.12987 


-0.14810 




LUMO 


-0.11542 


-0.12583 




gap/eV 


0.39 


0.61 




E 


-2952.70120 


-2952.70002 




AE 


46.7 


47.4 




Nimag 


none 


none 



The overall energy order for the (C4H6)2Fe2(CO) structures investigated in this section is 
lQ-1 < lT-1 < lQ-2 < lT-2 < lS-l-'lS-2 (M06-L). Thus the higher spin states are energetically 
favored. The energy gaps between the HOMOs and LUMOs in Tables 4, 5, and 6 for lS-1 lT-1 
lQ-1 are 0.39, 1.39 and 2.31 eV, respectively, and that for lS-2 lT-2 lQ-2 are 0.61, 1.93 and 
1.72 eV, respectively. 

3.3. NBO Analysis 

The natural charges on the iron atoms and the Wiberg Bond Indices (WBIs) for the iron-iron bonds 
are listed in Table 7 along with the Fe-Fe distances, the iron electronic configurations, and formal 
iron-iron bond orders. For the less unsaturated binuclear butadiene iron carbonyls 
(C4H6)2Fe2(CO)„ {n = 7, 6, 5, 4, 3) studied previously [4] the natural charges on the iron atoms 
correlate mainly with the numbers of carbonyl groups on the iron atoms and the WBIs correlated with 
the formal iron-iron bond orders. However, the much greater variety of iron electronic configurations 
and spin states encountered in the highly unsaturated (C4H6)2Fe2(CO)„ = 2, 1) structures reported in 
this paper make the interpretations much less clear. In this connection, iron atoms bonded to two CO 
groups were found to be essentially neutral. In most cases iron atoms bonded to only a single CO 
group (or half of two bridging carbonyl groups) have natural positive charges in the range 0.13 to 0.25. 
Carbonyl-free iron atoms are even more positive but their natural positive charges span a wide range 
from 0.28 to 0.83. Thus the previously observed [4] general trend of increasing natural negative 
charges on the iron atoms with increasing number of carbonyl ligands is also found here. However, 
other factors besides the number of carbonyl groups also affect significantly the natural charges on the 
iron atoms. 

The WBIs for the iron-iron bonds in the (C4H6)2Fe2(CO)„ (n = 2, 1) structures also showed the 
expected correlation of increased WBI with an increase in formal bond order (Table 7). However, the 
ranges of WBIs for a given formal iron-iron bond order are relatively broad indicating the significant 
influence of other factors. The one example of an Fe-Fe single bond in 2T-2 has WBI of 0.19. Most of 
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the Fe=Fe double bonds have WBIs in the range 0.26 to 0.43. However, there are some unusually high 
WBIs for apparent Fe=Fe double bonds including 0.59 for the doubly bridged coaxial (C4H6)2Fe2(CO)2 
structure 2S-2 and the very high value of 0.93 for the perpendicular (C4H6)2Fe2(CO)2 structure 2S-3. 
Interpreting the iron-iron bonds in these two structures as formal quadruple bonds would rationalize 
these significantly higher WBIs and give both iron atoms the favorable 18-electron configurations but 
would be inconsistent with the iron-iron distances of -23 A. The formal Fe=Fe triple bonds in the 
(C4H6)2Fe2(CO)^ {n = 2, 1) structures exhibit WBIs in the broad range 0.38 to 0.81. The two examples 
of formal Fe-Fe quadruple bonds, namely the Fe-Fe bonds in the singlet (C4H6)2Fe2(CO) structures 
lS-1 and lS-2, exhibit by far the highest WBIs at 1.30 and 1.42, respectively, consistent with the high 
formal bond orders. 

Table 7. Fe-Fe distances, NPA natural charges, iron electron configurations, traditional 
formal Fe-Fe bond orders and Wiberg bond indices (WBIs) for the (C4H6)2Fe2(CO)„ 
(n = 2, 1) structures using the BP86 method. Global minima structures are in bold type. 

Fe Natural Fe Electron Formal Wiberg 

Fe-Fe 

(11 -C4H6)2Fe2(CO)„ Charge Configuration Bond Bond 

Distance 











Fel 


Fe2 


Fel 


Fe2 


Order 


Index 


n = 2 


2Q-1 


Cs 


2.361 


0.164 


0.685 


17 


15 


2 


0.26 




2Q-2 


Ci 


2.252 


-0.066 


0.834 


17 


15 


2 


0.45 




2Q-3 




2.366 


0.467 


0.467 


16 


16 


2 


0.38 




2T-1 


Q 


2.197 


0.251 


0.251 


17 


17 


3 


0.52 




2T-2 


Ci 


2.441 


0.023 


0.335 


17 


17 


1 


0.19 




2T-3 




2.267 


0.157 


0.157 


17 


17 


3 


0.52 




2S-1 


Ci 


2.348 


0.095 


0.097 


18 


18 


2 


0.26 




2S-2 




2.332 


0.121 


0.121 


16 


16 


2 


0.59 




2S-3 


C2 


2.337 


0.206 


0.206 


16 


16 


2 


0.93 


n=l 


lQ-1 


Ci 


2.252 


0.187 


0.742 


17 


15 


3 


0.52 




lQ-2 


Ci 


2.262 


0.371 


0.663 


16 


16 


3 


0.38 




lT-1 


Cs 


2.183 


0.177 


0.682 


17 


15 


3 


0.86 




lT-2 


Cs 


2.108 


0.194 


0.503 


16 


16 


3 


0.81 




lS-1 


Cs 


2.051 


0.131 


0.467 


18 


16 


4 


1.30 




lS-2 


Cs 


2.041 


0.278 


0.385 


18 


16 


4 


1.42 



3.4. Vibrational Frequencies 

Table 8 exhibits the v(CO) frequencies and their infrared intensities for all of the (C4H6)2Fe2(CO)„ 
(n = 2, 1) structures, evaluated using the BP86 method, which has been shown to be a reliable 
predictor of such v(CO) frequencies. The terminal v(CO) frequencies fall in the range from 1908 to 
1962 cm"^ whereas the bridging v(CO) frequencies are significantly lower falling in the range from 
1761 to 1858 cm"^ (Table 8). The significantly lower v(CO) frequencies for bridging relative to 
terminal carbonyls is well-established and is consistent with the lower effective C-O bond order for 
bridging relative to terminal carbonyl groups in a given type of metal carbonyl structure. 
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Table 8. The v(CO) frequencies (cm"^) and their infrared intensities (km/mol, in 
parentheses) for (C4H6)2Fe2(CO)„ {n = 2, 1) structures as determined by the BP86 method. 
Bridging v(CO) frequencies are in bold type. 







1793^3 491^ 1942ra 1037^ 

A / \(*^ A ± VT^^V^Cl, ± / ^ 




2S-2 (Ci) 


1761 (bu, 883), 1783 (ag, 0) 




2S-3 (C2) 


1925 (b, 424), 1955 (a, 1040) 




2T-1 (C2) 


1842 (b, 1078), 1858 (a, 65) 




2T-2 (Ci) 


1908 (a, 485), 1957 (a, 1296) 




2T-3 (C2) 


1925 (a, 656), 1950 (a, 1043) 




2Q-1 (C) 


1815 (a', 454), 1962 (a', 1045) 




2Q-2 (Ci) 


1912 (a, 798), 1948 (a, 706) 




2Q-3 {C2h) 


1842 (b„, 1002), 1858 (ag, 0) 


(C4H6)2Fe2(CO) 


lS-1 (G) 


1933 (a', 901) 




lS-2 (Q) 


1826 (a', 716) 




lT-1 (Q 


1929 (a', 857) 




lT-2 (C) 


1835 (a', 681) 




lQ-1 (Ci) 


1925 (a, 870) 




lQ-2 (Ci) 


1802 (a, 666) 



3.5. Thermochemistry 

In order to check the potential experimental accessibility of the title compounds, we examined the 
following energies: 

(1) The dissociation energies of carbonyl groups from (C4H6)2Fe2(CO)2, namely: 

(C4H6)2Fe2(CO)2 ^ (C4H6)2Fe2(CO) + CO 

(2) The energies of the following disproportionation reaction: 

2(C4H6)2Fe2(CO)2 ^ (C4H6)2Fe2(CO)3 + (C4H6)2Fe2(CO) 

(3) The fragmentation energies of the binuclear (C4H6)2Fe2(CO)2 to mononuclear fragments by the 
following reaction: 

(C4H6)2Fe2(CO)2 ^ 2C4H6Fe(CO) 

Table 9 lists the energies and corresponding free energies for the above reactions taking the 
energies of the structures (C4H6)2Fe2(CO)3 and C4H6Fe(CO) from ref. [4]. The predicted energy for 
loss of a single carbonyl group from (C4H6)2Fe2(CO)2 is large, namely -30 kcal/mol as well as a free 
energy of --20 kcal/mol. The disproportionation of (C4H6)2Fe2(CO)2 to give (C4H6)2Fe2(CO)3 and 
(C4H6)2Fe2(CO) using the BP86 method is exothermic, while it is endothermic using M06-L and 
B3LYP methods by 3.4 and 15.6 kcal/mol, respectively. The corresponding free energies reveal a 
similar trend. The predicted fragmentation energies of (C4H6)2Fe2(CO)2 to mononuclear C4H6Fe(CO) 
is significantly larger, namely --60 kcal/mol by any of the three methods. 
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Table 9. Dissociation energy for removal of one carbonyl group from (C4H6)2Fe2(CO)2, 
disproportionation energy for (C4H6)2Fe2(CO)2, and dissociation energy for 

(C4H6)2Fe2(CO)2 2C4H6Fe(CO). The corresponding free energies are in italics (kcal/mol). 



BP86 M06-L B3LYP 



(C4H6)2Fe2(CO)2(2Q-l) ^ (C4H6)2Fe2(CO)(lQ-l) +CO 



2(C4H6)2Fe2(CO)2(2Q-l) ^ (C4H6)2Fe2(CO)3(3S-l) + 
(C4H6)2Fe2(CO) (lQ-1) 



(C4H6)2Fe2(CO)2(2Q-l) ^ 2C4H6Fe(CO) 



G 
E 
G 



G 



E 



E 



34.1 33.9 30.0 

22.1 22.5 15.6 

-13.9 3.4 15.6 

-10.5 9.2 18.0 

59.1 73.7 59.3 

46.6 62.6 47.6 



4. Conclusions 

Unsaturation in binuclear metal carbonyl derivatives can lead to metal-metal multiple bonding, 
four-electron donor bridging carbonyl groups, and/or metal electronic configurations less than the 
favorable 18-electron configurations. None of the highly unsaturated (C4H6)2Fe2(CO)„ (n = 2, 1) 
structures found in this work has a four-electron donor bridging carbonyl group. Instead the lowest 
energy (C4H6)2Fe2(CO)„ (n = 2, 1) structures are perpendicular structures having iron atoms 
with 15- and 17-electron configurations. This leads to quintet spin states in addition to iron-iron 
multiple bonds of formal order two for (C4H6)2Fe2(CO)2 and three for (C4H6)2Fe2(CO). In addition, 
agostic hydrogen atoms forming C-H-Fe bridges are seen to be a feature of (C4H6)2Fe2(CO)„ (n = 2, 1) 
structures, albeit not the lowest energy such structures. In such structures a butadiene ligand is bonded 
to one of the iron atoms as a tetrahapto ligand and to the other iron atom through two C-H-Fe units 
from the end CH2 groups with agostic hydrogen atoms bridging iron-carbon bonds. Singlet 

o 

(C4H6)2Fe2(CO) structures with formal Fe-Fe quadruple bonds of lengths --2.05 A were also found but 
at very high energies (--47 kcal/mol) relative to the global minimum. 
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